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The invention concerns a simulation method for 
time dependent forecast of a detailed composition with 
Q components and/or pseudo-components of a fluid 
produced by a reservoir, which consists in: 
representing the reservoir in the form of a grid system 
(j) whereof each grid constitutes an elementary volume 
filled with fluid; producing, in a manner known per se, 
a compositional simulation of the aggregate fluid with 
N components and pseudo-components (i), N being less 
than Q. The method is characterised it that it further 
consists in: e) determining, at each time step (m) and 
for each grid (j) tne liquid and fluid vapour phase 
composition for detailed modelling with Q components 
and pseudo-components, on the basis of the vaporised 
fraction (.theta.jm) and the equilibrium constants 
(kijm) of the aggregate fluid; f) evaluating for each 
grid, of time step (m+1), the molar amount of each of 
the Q components and/or pseudo-components of the 
itemised fluid on the basis of the values corresponding 
to time step (m), of the flow rates of the aggregate 
fluid phases and the liquid and itemised fluid vapour 
phase compositions determined at step e); and g) 
evaluating, for each production well, the detailed 
composition of the fluid produced between times t and 
t' corresponding to time steps m and m+1, from the 
aggregate fluid flow rates and the liquid and itemised 
fluid vapour phase compositions determined at step e). 
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(57) Methode de simulation pour predire en fonction du 
temps une composition detail lee a Q compos ants et/ou 
pscudo-eomposanls d'un fluidc produiL par un reservoir. 
Elle est du type consistant a: representor le reservoir sous 
la forme d'un reseau de mailles (j) dont chacune 
constitue un volume elementaire rempli de fluide, 
realiser, de facon connue en soi, une simulation 
oornposilionnelle du fluide regroupe a N eomposants cL 
pseudo- eomposants (i), N etant iiiferieur a Q, et elle est 
caracterisee en ce qu'elle consiste en outre a: e) 
determiner, a chaque pas le temps (m) et pour chaque 
maille (j) la composition des phases liquide et vapeur du 
fluidc pour la modelisation dclaillec a Q eomposanls el 
pseudo-eomposants, a partir des valeurs de la fraction 

vaporisee (9j m ) et des constantes d'equilibre (kj 111 ) du 

fluide regroupe, f) evaluer pour chaque maille, au pas de 
temps (m+1), la quantite molaire de cliacun des Q 
compos ants et/ou pseud o-composants du fluide detaille a 
partir des valeurs correspondantes au pas de temps (m), 
des debits des phases du fluide regroupe et des 
compositions des phases liquide et vapeur du fluide 
detaille determinees a Tetape e), et a g) evaluer, pour 
chaque putts de production, la composition detaillee du 
fluide produit entre des instants t et t' correspondant aux 
pas de temps m et m+1 3 a partir des debits des phases du 
fluide regroupe et des compositions des phases liquide el 
vapeur du fluide detaille determinees a l'etape e). 



(57) The invention concerns a simulation method for time 
dependent forecast of a detailed composition with Q 
components and/or pseudo-components of a fluid 
produced by a reservoir, which consists in: representing 
the reservoir in the form of a grid system (j) whereof each 
grid constitutes an elementary volume filled with fluid; 
producing, in a manner known per se, a compositional 
simulation of Lhe aggregate fluid with N components and 
pseudo-components (i), N being less than Q. The method 
is characterised it that it further consists in: e) 
determining, at each time step (m) and for each grid (j) 
the liquid and fluid vapour phase composition for 
detailed modelling with Q components and pseudo- 
components, on the basis of the vaporised fraction (9- m ) 

and the equilibrium constants (k-- m ) of the aggregate 

fluid; f) evaluating for each grid, of time step (m+1), the 
molar amount of each of the Q components and/or 
pseudo-components of the itemised fluid on the basis of 
the values corresponding to time step (m), of the flow 
rates of the aggregate fluid phases and the liquid and 
itemised fluid vapour phase compositions determined at 
step e); and g) evaluating, for each production well, the 
detailed composition of the fluid produced between 
times t and t' corresponding to time steps m and m+1, 
from the aggregate fluid flow rates and the liquid and 
itemised fluid vapour phase compositions determined at 
step e). 
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(57) Abstract 



The invention concerns a simulation method for time dependent forecast of a detailed composition with Q components and/or 
pseudo-components of a fluid produced by a reservoir, which consists in: representing the reservoir in the form of a grid system (j) whereof 
each grid constitutes an elementary volume filled with fluid; producing, in a manner known per sc, a compositional simulation of the 
aggregate fluid with N components and pseudo-components (i), N being less than Q. The method is characterised it mat it further consists 
in: e) determining, at each time step (m) and for each grid (j) the liquid and fluid vapour phase composition for detailed modelling with Q 
components and pseudo-components, on the basis of the vaporised fraction (0/™) and the equfflbrium constants (ky^) of the aggregate fluid- 
f) evaluating for each grid, of time step (m+1), the molar amount of each of the Q components and/or pseudo-components of the itemised 
fluid on the basis of the values corresponding to time step (m), of the flow rates of the aggregate fluid phases and the liquid and itemised 
fluid vapour phase compositions determined at step e); and g) evaluating, for each production well, the detailed composition of the fluid 
produced between times t and t corresponding to time steps m and m+1, from the aggregate fluid flow rates and the liquid and itemised 
fluid vapour phase compositions determined at step e). 

(57) Abreg* 



Methode de simulation pour pnSdire en fonction du temps une composition d&aillee a Q composants et/ou pseudo-composants d'un 
fluide produit par un reservoir. Elle est du type consistant a: representer le reservoir sous la forme d'un reseau de mailles (j) dont chacune 
constitue un volume 616mentaire rempli de fluide, realiser, de racon connue en soi, une simulation compositionnelle du fluide regroupe" a 
N composants et pseudo-composants (i), N etant infeneur a Q, et elle est characterised en ce quelle consist© en outre a: e) determiner, a 
chaque pas Ie temps (m) et pour chaque maille (j) la composition des phases liquide et vapeur du fluide pour la mode*lisation d6taill6e a*Q 
composants et pseudo-composants, a partir des valeurs de la fraction vaporised (0] m ) et des constantes d'equilibre (kj m ) du fluide regroupe' 
f) evaluer pour chaque maille, au pas de temps (m+l) f la quantity molaire de chacun des Q composants et/ou pseudo-composants du fluide 
d6taill6 a partir des valeurs correspondantes au pas de temps (m), des debits des phases du fluide regroupe et des compositions des phases 
liquide et vapeur du fluide dtftaille* determiners a l*6tape e), et a g) 6valuer, pour chaque puits de production, la composition d6taillee du 
fluide produit entre des instants t et t' correspondant aux pas de temps m et m+l, a partir des dibits des phases du fluide regroup* et des 
compositions des phases liquide et vapeur du fluide ddtaille" dererminees a I'etape e). 
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Simulation method for time dependent forecast of 
detailed composition of a fluid produced by a reservoir 

The present invention relates to a method of simulation 
5 for predicting, as a function of time, a detailed 
composition of a fluid produced by a reservoir and more 
particularly a detailed composition of a fluid 
contained in and produced by an oilfield in which one 
or more production wells are sunk. 

10 

The compositional simulation of an oilfield is commonly 
used to provide forecast production profiles for the 
field which make it possible in particular to determine 
the most suitable production scheme for this field. 

15 

The compositional simulation of an oilfield is 
implemented not by using a real description of the 
fluid but by using a fluid modelled by a number of 
components which is smaller than the number of 

20 components of the real fluid. Specifically, the number 
of components of the reservoir fluid being relatively 
large, a model using all the components would lead to 
excessive calculation times. This \ prohibitive 
calculation time has prompted specialists to group the 

25 pure components into pseudo-components , for example a 
pseudo-component which groups together nitrogen and 
methane, a pseudo-component which groups together C 3 
and C 4 hydrocarbons, etc., and to carry out the 
simulation of the field on a reduced number N of pure 

30 components and of pseudo-components. Such a reduced 
composition is termed a "grouped composition of the 
fluid" or "grouped fluid". Generally N represents 5 to 
10 components and pseudo-components, this being 
regarded as sufficient for a proper representation of 

35 the behaviour of the reservoir fluid under bottom 
conditions. The grouped fluid is described by one of 
the equations of state which are well known to 
specialists and is adapted to the said grouped fluid. 
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Before performing the compositional simulation of the 
field, the reservoir or field is represented in the 
form of a grid of ceils each of which constitutes an 
elementary volume of the said reservoir. The number of 
5 cells may be as many as several thousand and each cell 
exhibits properties specific to it, such as geometry, 
porosity, permeability, etc. Moreover, at least one 
injection well or production well which is sunk into 
the said reservoir may correspond to certain cells. 

10 

Simulation of the field makes it possible to calculate 
for each cell a number of principal variables of the 
reservoir fluid, the said variables possibly being the 
quantity (number of moles) of fluid, the composition or 

15 molar fraction of each component and pseudo-component 
of the grouped fluid and the pressure prevailing in 
each cell. These variables are known at the instant t=C 
(start of the exploiting of the field) and are 
calculated thereafter by the field simulation at each 

20 instant t or timestep m. For each cell, it is possible, 
on the basis of these principal variables, to calculate 
any other property of the fluid present in the cell, 
such as the number of phases, the composition of each 
phase, etc. For each cell associated with a production 

25 well, it is possible, on the basis of these principal 
variables and of the production constraints imposed, to 
also calculate the production flow rate and the 
composition of the fluid produced by the said well. 

30 The thermodynamic properties of the grouped fluid can 
be calculated using the equation of state associated 
with it. 

Modelling with N components, although it is sufficient 
35 to represent the behaviour of the reservoir fluid under 
bottom conditions, is no longer appropriate for 
simulating the behaviour of the fluid in plants for its 
surface exploitation, the latter requiring knowledge of 
a more detailed composition with Q components and/or 
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pseudo-components of the fluid produced by the field, Q 
being greater than N and for example of the order of 16 
to 30 . 

5 Hitherto, this detailed composition with Q components 
and/or pseudo-components was obtained on the basis of 
the grouped composition N by assuming that the 
composition of the pseudo-components remained constant 
over time. 

10 

Such a procedure gives rise to errors since the 
composition of the pseudo-components varies over time, 
especially as a function of the pressure of the field 
or following the injection of a gas into the field. 

15 

Therefore, the detailed composition obtained remains 
very approximate and does not allow satisfactory 
prediction of the behaviour of the fluid under the 
conditions of surface exploitation several years after 
20 the start of production. 

The purpose of the present invention is to propose a 
method of simulation predicting a detailed composition 
of the fluid produced by a reservoir while taking into 

25 account the parameters of the reservoir and which is 
much more accurate, such a prediction being dubbed 
"delumping" by specialists. v De lumping" consists in 
predicting what the results of a reservoir simulation 
would have been if the detailed modelling of the fluid 

30 had been used, and to do so on the basis of the results 
of a simulation carried out by using the grouped 
modelling of the said fluid. 

The subject of the invention is a method of the type 
35 consisting in: 

a) representing the reservoir in the form of a grid of 
cells (j) each of which constitutes an elementary 
fluid- filled volume , 
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b) defining the fluid by a grouped modelling with N 
components and pseudo-components (i) and determining an 
equation of state describing the fluid in this grouped 

5 modelling, 

c) also defining the fluid by a detailed modelling with 
Q components and/or pseudo-components, Q being greater 
than N, and determining another equation of state 

10 describing the fluid in this detailed modelling, 

d) carrying out, in a manner known per se, a 
compositional simulation of the grouped fluid with N 
components and pseudo-components (i) , the said 

15 compositional simulation making it possible to 
calculate at least for each cell (j) and at consecutive 

timesteps (m, m+1, ...) the vaporized fraction (©j), 

m 

the liquid/vapour equilibrium constants ( k i j ) for each 
component (i), the injection or production flow rates 

m 

20 ( s i j ) and for each pair of cells (j, h) the flow rates 

Etl m 

of the liquid phase (uojh) and vapour phase (ugjh) of 
the fluid with N components and pseudo-components, and 
it is characterized in that it furthermore consists in: 

25 e) determining, at each timestep (m) and for each cell 
(j) the composition of the liquid and vapour phases of 
the fluid for the detailed modelling with Q components 
and/or pseudo-components, on the basis of the values of 

the vaporized fraction (6 j ) and of the equilibrium 

m 

30 constants ( k i j ) of the grouped fluid, 

f ) evaluating for each cell, at timestep (m+1) , the 
molar quantity of each of the Q components and/or 
pseudo-components of the detailed fluid on the basis of 
35 the corresponding values at timestep (m) , of the flow 
rates of the phases of the grouped fluid and of the 
compositions of the liquid and vapour phases of the 
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detailed fluid which were determined in stage e) , and 
in 



g) evaluating, for each production well, the detailed 
5 composition of the fluid produced between instants t 
and t' corresponding to timesteps m and m+1, on the 
basis of the flow rates of the phases of the grouped 
fluid and of the compositions of the liquid and vapour 
phases of the detailed fluid which were determined in 
10 s~age e) . 



According to another characteristic of the present 
invention, stages e, f and g are implemented at the 
same time as the carrying out of the compositional 
15 simulation of the grouped fluid. 

According to another characteristic of the present 
invention, the results of stage d are stored in memory 
and then used later for the implementation of stage e, 
20 f and g. 



The method according to the present invention is 
implemented in respect of a reservoir consisting, for 
example, of an oilfield. A field simulation is carried 
25 out in a known manner. To this end, the field is 
represented in the form of a grid of cells, some of the 
cells or a group of cells being associated with 
production wells sunk into the oilfield to be 
exploited . 

30 

To reduce the calculation time, the compositional 
simulation of the field is performed on a limited 
number N of components and pseudo-components, for 
example from 5 to 10, defined in the manner indicated 
35 earlier, these components and pseudo-components being 
selected as a function of the nature of the field. 
As indicated earlier, the grouped fluid with N 
components and pseudo-components is described by an 
equation of state which could for example be the 
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PENG-ROBINSON equation adapted to the said grouped 
composition . 

The compositional simulation of the field makes it 

5 possible to calculate, at each timestep and for each 

cell, a number of principal variables of the grouped 

fluid, which variables, in this case, are: 

m 

the quantity of fluid: f j 

m 

the molar fraction of component i: z i j with 1 < i < N 

m 

10 the pressure: pj 



Through prior measurements, the value of the said 
principal variables is known at the start of the 
exploiting of the field, that is to say at the instant 
15 t=0 or at timestep m=0 . 

To go from timestep m to timestep m+1, the 
compositional simulation comprises several stages which 
are detailed hereinbelow. 

20 

In a first stage and for timestep m, a liquid/vapour 
equilibrium (flash) calculation is performed on the 
grouped fluid, for each cell j so as to determine: 

A _ . ... _ .. j 

the molar fraction of each component i in the liquid 

ru 

phase (xi with i = 1...N) and 

30 the molar fraction of the said component i in the 

m 

vapour or gas phase <yi with i = 1...N). 



These various fractions are determined with the 
following system of equations: 

35 

°i(*,7 - ^(y™ ygj) « isi* n 
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5 



25 



30 



(i-ejOxjjUejvy-ij 1 



and 

N N 



in which: 

<T>i is the fugacity of component i 



The vaporized fraction 9™ and the equilibrium constants 
m 



10 k™ = — jjp with i varying between 1 and N are stored in 
x ij 

memory so as to use them in another stage dubbed 
"delumping" which will be explained later. 

In a second stage and for each cell j, certain of the 
15 properties of each liquid and gas phase are evaluated, 

.on 



namely the density (p^ and Pgj), viscosity, saturatic 



and relative permeability. 

In a third stage, for each cell j and at timestep m, we 
20 calculate the coefficients of the equation for the 
pressure . 

The pressure equations for a cell j at the instant 
t + At or at timestep m + 1 are: 



At 



- V T m fn m+1 -n ra+1 ) + ? m 

bfi-cD Jh h j j j 



To solve these equations, the coefficients to be 
calculated are: 

- the total compressibility C™ (fluid + rock) , 

- the generalized transmissivities 
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TS = Z«y l (^ i +X g y i ) ra t jh 

i=l 



- the source term: 

N 

In these equations: 

h are the cells neighbouring the cell j, h£ J{j), 

a^. is the partial molar volume of component i in cell 

j at timestep m, X,™j and A,^ are the mobilities of the 

10 liquid (o) and gas (g) phases, t jh is the transmissivity 
between the centres of cells j and h, 1 being the 
upstream cell. 

In a fourth stage, the pressure equations are solved 
15 f or all the cells at the instant t + At, that is to say 
at timestep m + 1 . 

These equations form a linear system of the type: 

2 0 A.p=B 

« « « 

f \ T 

where p= |p^ +1 ,...pj +1 j , J being the number of cells of 
the grid. 

These equations are solved by a standard iterative 

2 5 procedure. 

In a fifth stage, the flow rates of the phases between 
each pair of adjacent cells and for each associated 
well are calculated at timestep m + 1 on the basis of, 

3 0 for example, DARCY' s law. 



CA 02286173 199910-06 



- 9 - 



u 



m 
ojh 



gjh 



rn 



Poj x ol- l jhVPh Pj ) 



p being the density and 1 being the upstream cell with 
1 = j or h. 



According to one characteristic of the invention, the 
flow rates of the phases which will be used in the 
delumping are stored or kept in memory. 

10 In a sixth stage, the quantity f and the composition z 
of the grouped fluid are evaluated, in each cell j and 
at timestep m + 1 , on the basis of the values at 
timestep m, of the flow rates of the phases which were 
evaluated in the fifth stage and of the compositions of 

15 the phases which were calculated in rhe first stage. To 
do this, the conservation laws defined by the following 
equations are applied: 



20 Through an iterative process, the above stages are 
repeated so as to calculate the principal variables at 
times t + At, t + 2 At, t + 3At, etc. or else at the 
corresponding timesteps m + 1, m + 2, m+3 etc, the 
increment At being constant or variable. 





Since the flow can occur in any direction between 
neighbouring cells, all the cells must be processed at 
each timestep. 
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According to the present invention, a detailed 
modelling with Q components and/or pseudo-components is 
defined together with another equation of state 
describing the fluid thus modelled. 

5 

To do this, the components and/or pseudo-components 
whose characteristics one wishes to know and which are 
necessary in order to represent the production fluid 
under the surface conditions are chosen beforehand. The 
10 number Q which is greater than the number N of the 
grouped composition is generally of the order of 16 to 
30. The composition of the detailed fluid with Q 
components and/or pseudo-components is known at 
timestep m = 0 . 

15 

The data stored in the first and fifth stages are used 
to extend them to the Q components and/or 
pseudo-components of the detailed composition, in such 
a way as to perform the " delurrtping" operation. For the 
20 '"delumping" operation it is assumed that the reservoir 
fluid comprises just two non-aqueous phases (liquid and 
gas) and that the water plays no role in the liquid/gas 
equilibrium and that none of the other components of 
the fluid dissolves in the water. 

25 

In what follows, upper case letters are used to 
represent the detailed fluid and lower case letters are 
used to represent the grouped fluid, 

30 To carry out the ,% de lumping" , the further assumption is 
made that at each timestep m, 

the molar fraction of the vaporized phase in each 
cell is independent of the number of constituents 
chosen to model the fluid, 

35 




the number of moles of liquid and of vapour 
flowing between the cells (U, u) and into or out 
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of the wells (S, s) is independent of the number 
of constituents chosen to model the fluid. 



10 



To perform the " delumping" , it is necessary to know the 
equilibrium constants K for the detailed fluid as well 
as the molar fraction © already determined through the 
equality 0=9. 

The calculation of the equilibrium constants can be 
performed on the basis of the general equation: 



Ln[K I ] = AC 0+ t AC n .ir in ou AC n -C^-Cjf 

15 

as indicated by C 
Consistent Procedure 
published in FLUID 
225-232. 

20 

For a two-parameter equation of state such as the PENG- 
ROBINSON equation, the generalized LEIBOVICI equation 
can be written for the detailed fluid in the following 
manner, so long as there is no binary interaction 
25 coefficient: 



LEIBOVICI in his article "A 
for Pseudo-Component Delumping'' 
PHASE EQUILIBRIA, 117 (1996), 



where ai and hj are the parameters of the equation of 
state for the component I and where the parameters 
30 rj ™ . j3 ^ and y^ 1 are parameters calculated from the 



equilibrium constants 




for the grouped fluid by 



minimizing the function: 
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f(if .Pf ,rf ) = fjn? A + rf b, - u(k»)J 



The parameters 7 J^Jand^ thus obtained will be 

approximately the same for both modellings, detailed 
and grouped, of the fluid, the said parameters 
depending only on the phase parameters, the pressure 
and the temperature. 



In the presence of a binary interaction coefficient, 
the parameters t|, p and y can be determined by 
regression . 



Therefore, when the detailed composition is known 

in cell j at timestep m, the detailed compositions of 
15 the liquid and vapour phases can be estimated in the 
same cell j and for the same timestep m. 

I being the component and/or pseudo-component of the 
detailed composition. These compositions can be 
20 normalized so that their sum is equal to 1 if 
necessary . 

When the compositions of the detailed phases are known 
for all the cells of the grid at timestep m, it is 
25 possible to estimate the detailed compositions at 

the next timestep m + 1 from the equations below: 



ZffFF-^YjpsS + XfSj)-* V (y^U^X^) 
H _ — ^1 with I-l Q 



where 



Ff +l = Tf- At(s5 +SJJ) - At T (ug h + ug h ) 

helfj) 



30 in which j 1 = j for a flow rate from ceil j to cell h 
or into the well and j T = h for a flow rate from cell h 
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to cell j, and j' corresponding to the fluid injected 
in the case of injection wells, S then being negative, 
and F is the molar quantity of fluid in each cell with 
F = f . 

The molar flow rate for detailed component 1 produced 
by well W at timestep m is given by: 



10 and the detailed composition of the fluid produced by 
the well is obtained by normalizing the above flow 
rates : 



15 

By virtue of the present invention, a detailed 
composition with Q components and/or pseudo-components 
is obtained which enables the forecast production 
profiles to be better defined with enhanced accuracy as 

20 compared with the prior methods* This greater accuracy 
of the forecast production profiles allows a more 
reliable choice of development scheme (type and size of 
surface plants, number of wells, etc.), thus resulting 
in an increase in profitability and a decrease in 

25 economic risk. 



In the foregoing, certain values obtained in the course 
of the field simulation are stored or held in memory so 
as to be used later in the delumping. It goes without 
30 saying that these same values could be used directly if 
it were decided to perform the delumping at the same 
time as the said field simulation. 



Another advantage of the present invention is that the 
35 method can be used to calculate the concentration of 
hydrocarbons present in trace amounts in the fluid 
produced without having to perform a reservoir 
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simulation incorporating these hydrocarbons, provided 
that their equation of state parameters are known. 

It should foe noted that it is possible to treat the 
5 cases where the initial composition of the fluid varies 
with depth or laterally inside the reservoir, on 
condition that the initial detailed composition is 
known at every point of the reservoir. 

10 Two examples are given below to demonstrate the speed 
and accuracy of the method according to the invention, 
the tests having been performed on a parallepipedal 
reservoir model comprising four layers of various 
permeabilities, with a dip in the x direction, 

15 

The predictions were calculated so as to determine the 
behaviour of the fluid under the surface conditions 
over a period of 30 years from the start of 
exploitation, the fluid consisting of a volatile oil. 

20 

In the first example, the reservoir includes a 
water /oil contact, whilst in the second example the 
reservoir is devoid of water and contains only oil and 
gas . 

25 

In the modelling of the detailed fluid, the fluid is 
defined by 16 components and/or pseudo-components, as 
indicated in Table 1, and the binary interaction 
coefficients are regarded as equal to zero. 

30 

In the modelling of the grouped fluid, the fluid is 
defined by 7 pseudo-components, as indicated in 
Table 2. 

35 A- the temperature of the reservoir, the bubble-point 
pressure of the volatile oil is 338.5 bar for the 
detailed fluid and 336.6 bar for the grouped fluid- 
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Example 1 

The initial pressure of the reservoir is 400 bar and 
the reservoir is depleted through a single well which 
is drilled near the top of the reservoir. The oil 
5 production from the well is 300 m 3 per day and the 
maximum production of gas is 10 6 m 3 normal per day. 

A first compositional simulation of the reservoir was 
initially performed using the detailed modelling of the 

10 fluid with 16 components and/or pseudo-components with 
the aid of the equations of state adapted to each of 
the components. At the start, the well produces oil at 
the abovementioned rate. After around 1250 days, the 
bubble-point pressure is reached. The flow rate of gas 

15 increases continually up to the maximum value of 
10 6 N m 3 /day which is reached after around 5000 days. 
The well discharges at this maximum whilst the flow 
rate of oil decreases until the bottomwell pressure 
limit is reached at around the 6500 th day. After which, 

20 the two flow rates, oil and gas, decrease, the 
simulation terminating at the end of 10950 days (30 
years) . At this juncture the average pressure of the 
reservoir is slightly above 50 bar. Water production 
occurs after 5000 days. 

25 

A second compositional simulation of the reservoir is 
then repeated with the modelling of the grouped fluid 
with 7 pseudo-components. The production profiles 
achieved are substantially the same as those of the 
30 simulation with the detailed fluid. This suggests that 
the modelling with 7 pseudo-components is adequate to 
model the phase behaviour of the fluid under the 
conditions of the reservoir, at the very least for the 
purposes of studying the field. 

35 

When "delumping" is applied to the results obtained 
with the second simulation performed on the grouped 
fluid, it is noted that there is agreement between the 
results obtained after "deluraping" according to the 
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invention and those obtained with the first simulation 
which is performed with the aid of the equations of 
state adapted to each of the 16 components and/or 
pseudo-components. The difference lies in the fact that 
5 the first simulation of the detailed fluid requires a 
calculation time of 730 seconds whilst the overall 
calculation time for performing the second simulation 
(232 s) and the delumping (27 s) is only 259 seconds, 
that is to say the calculation time is reduced by 65% 
10 while still having excellent results. 

Example 2 

In this example, at the start the reservoir is 
15 completely filled with oil. The oil is produced by a 
well opening into the lower part of the reservoir, with 
a flow rate of 1000 m 3 /day. 

Gas is injected through a well, into the upper part of 
20 the reservoir, the said injection gas not consisting of 
a mixture of grouped components, but comprising, by 
molar percentage, 80% C lf 15% C 2 and 5% C 3 , this 
composition being used for the simulation with the 
detailed fluid and also for the "delumping". For the 
25 reservoir simulation with the grouped fluid, the 
composition of the injected gas is represented in molar 
percentages by 80% N 2 Ci, 15% C0 2 C 2 and 5% C 3 C 4 ; the 
injection of gas is carried out at a rate of 500 m 3 /day 
for the first 15 years and 900 m 3 /day for the next 
30 15 years. 

As in Example 1, a first simulation with a detailed 
fluid model with 16 components is performed* At the 
start, the flow rate of the injection gas is 
35 insufficient to maintain the pressure of the reservoir; 
therefore, the bubble-point pressure is quickly reached 
and the gas/oil ratio increases. After 15 years, when 
the pressure of the reservoir is 120 bar, the flow rate 
of the injection gas is increased so that after 



CA 02286173 199910-06 



- 17 - 

30 years, at the end of the simulation, the pressure of 
the reservoir is around 100 bar. 

A second simulation with the model of the grouped fluid 
is performed. 

The same effects as in Example 1 are noted, that is to 
say the model with 7 pseudo-components is still 
adequate to model the phase behaviour of the fluid 
under the conditions of the reservoir. 

When "delumping" according to the invention is applied 
to the results of the second simulation, it is noted 
that, by virtue of the invention, the calculation time 
for performing the second simulation (173 s) and the 
"delumping" (18 s) is 191 seconds, this representing 
around 15% of the calculation time required to perform 
the first simulation with the model of the detailed 
fluid, the results being excellent, as in Example 1. 

When binary interaction coefficients are introduced, it 
is noted that the agreement is still satisfactory, even 
though it is not of the same order as that of 
Examples 1 and 2. For the heaviest components, the 
agreement is almost complete. Only the C0 2 ~"C 2 
pseudo-components of the grouped fluid gave relatively 
weak agreement. On the other hand, the gain with regard 
to the calculation time is still very appreciable since 
it represents only 35% of the calculation time of the 
simulation performed on the detailed fluid model with 
equations of state. 
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TABLE 1 



Component 


Initial oil coraDosition finnl fr^n \ 


N 2 


0 00080 


C0 2 




Ci 


0.64 641 


c 2 


0 04165 


c 3 


0.04649 


IC 4 


0.01569 


NC 4 


0 . 02690 




0 . 01024 


NC 5 


0.00948 


c 8 


0 . 01915 


c 7 


0. 02003 


c 8 


0.01925 


c 9 


0.01391 




0. 01120 


CNGi* 


0. 06644 


CNH 2 * 


0. 02309 



TABLE 2 

5 



Component 


Initial oil composition (mol . frac.) 


N 2 Ci 


0. 64721 


C0 2 C 2 


0. 07094 


C3C4 


0. 08908 


C5C7 


0. 05890 


C 8 C ao 


0. 04435 


CNG X . 


0 . 06644 


CNH 2 - 


0.02309 



*CNGi and CNH 2 represent the heavy part of the fluid 
(Cn + ) , which part is divided into two 

pseudo-components, one (CNGi) found principally in the 
10 gas and the other (CNH 2 ) in the oil- 
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CLAIMS 

Method of simulation for predicting, as a function 
of time, a detailed composition of a fluid 
produced by a reservoir, of the type consisting 
in : 

a) representing the reservoir in the form of a 
grid of cells (j) each of which constitutes an 
elementary fluid-filled volume, 

b) defining the fluid by a grouped modelling with 
N components and pseudo-components (i) and 
determining an equation of state describing the 
fluid in this grouped modelling, 

c) also defining the fluid by a detailed modelling 
with Q components and/or pseudo-components, Q 
being greater than N, and determining another 
equation of state describing the fluid in this 
detailed modelling, 

d) carrying out, in a manner known per se, a 
compositional simulation of the grouped fluid with 
N components and pseudo-components (i), the said 
compositional simulation making it possible to 
calculate for each cell (j) and at consecutive 
timesteps <m, m+1, , at least the vaporized 

fraction (9 j ) r the liquid/vapour equilibrium 

TCI 

constants ( K i j ) for each component (i) , the 

injection or production flow rates (sj) and for 
each pair of cells (j, h) the flow rates of the 
liquid phase (uojh) and vapour phase (ugjh) of 
the fluid with N components and pseudo-components, 
and it is characterized in that it furthermore 
consists in: 

e) determining, at each timestep (m) and for each 
cell (j) the composition of the liquid and vapour 
phases of the fluid for the detailed modelling 
with Q components and pseudo-components, on the 
basis of the values of the vaporized fraction 
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m m 

{U j ) and of the equilibrium constants ( k i j ) of 
the grouped fluid, 

f ) evaluating for each cell, at timestep (m+1) , 
the molar quantity of each of the Q components 
and/or pseudo-components of the detailed fluid on 
the basis of the corresponding values at timestep 
(rn) , of the flow rates of the phases of the 
grouped fluid and of the compositions of the 
liquid and vapour phases of the detailed fluid 
which were determined in stage e) , and in 

g) evaluating, for each production well, the 
detailed composition of the fluid produced between 
instants t and t' corresponding to timesteps m and 
m+1, on the basis of the flow rates of the phases 
of the grouped fluid and of the compositions of 
the liquid and vapour phases of the detailed fluid 
which were determined in stage e) . 

Method according to Claim 1, characterized in chat 
stages e, f and g are implemented at the same time 
as the carrying out of the compositional 
simulation of the grouped fluid. 

Method according to Claim 1, characterized in that 
the results of stage d are stored in memory and 
then used later for the implementation of stage e, 
f and g. 



